Maximum elastic deformations of relativistic stars 
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We present a method for calculating the maximum elastic quadrupolar deformations of relativistic 
stars, generalizing the previous Newtonian, Cowling approximation integral given by [G. Ushomirsky 
et al, Mon. Not. R. Astron. Soc. 319, 902 (2000)]. (We also present a method for Newtonian 
gravity with no Cowling approximation.) We apply these methods to the m — 2 quadrupoles most 
relevant for gravitational radiation in three cases: crustal deformations, deformations of crystalline 
cores of hadron-quark hybrid stars, and deformations of entirely crystalline color superconducting 
quark stars. In all cases, we find suppressions of the quadrupole due to relativity compared to the 
Newtonian Cowling approximation, particularly for compact stars. For the crust these suppressions 
are up to a factor ~ 6, for hybrid stars they are up to ~ 4, and for solid quark stars they are at 
most ~ 2, with slight enhancements instead for low mass stars. We also explore ranges of masses 
and equations of state more than in previous work, and find that for some parameters the maximum 
quadrupoles can still be very large. Even with the relativistic suppressions, we find that 1.4M@ 
stars can sustain crustal quadrupoles of a few x 10 39 g cm 2 for the SLy equation of state or close 
to 10 40 g cm 2 for equations of state that produce less compact stars. Solid quark stars of 1.4 Mq 
can sustain quadrupoles of around 10 44 g cm 2 . Hybrid stars typically do not have solid cores at 
1AM®, but the most massive ones (~ 2Mq) can sustain quadrupoles of a few x 10 41 g cm 2 for 
typical microphysical parameters and a few x 10 42 g cm 2 for extreme ones. All of these quadrupoles 
assume a breaking strain of 10 _1 and can be divided by 10 45 g cm 2 to yield the fiducial "ellipticities" 
quoted elsewhere. 



PACS numbers: 04.30.Db, 04.40.Dg, 97.60.Jd 



I. INTRODUCTION 

Shortly after the discovery of pulsars and the realiza- 
tion that they are rotating neutron stars, deformations of 
rotating neutron stars were proposed as sources of con- 
tinuous gravitational radiation jll-Q ; see [f| for an early 
review. Searches for such radiation are an ongoing con- 
cern of the LIGO and Virgo gravitational wave detec- 
tors see [§-[ll| for recent reviews. It is thus of great 
interest to know the maximum quadrupolar deformation 
that a neutron star could sustain, in order to motivate 
further searches and help interpret upper limits or de- 
tections. In the case of elastic (as opposed to magnetic) 
deformations, the main factor influencing the answer is 
whether the neutron star contains particles more exotic 
than neutrons 0, [l2j. However, the structure of the star 
also plays an important role. 

While there are relativistic calculations of the 
quadrupole deformations due to magnetic fields 
(e.g., [Hil - [l6l ]). all the computations involving elastic 
deformations have used Newtonian gravity. Moreover, 
all but two of these computations have used the integral 
expression obtained in the Cowling approximation 
(i.e., neglecting the self-gravity of the perturbation) 
by Ushomirsky, Cutler, and Bildsten (UCB) fl7j ; see 
[llElS]- Haskell, Jones, and Andersson (HJA) HJ 
dropped the Cowling approximation using a somewhat 
different formalism than UCB's; there is a further 
application of their results in (22j . 

We improve these treatments by generalizing the UCB 



integral to relativistic gravity with no Cowling approx- 
imation. We also provide a similar generalization for 
the Newtonian no-Cowling case, as a warm-up. In ad- 
dition to providing a simpler formalism for performing 
computations than the more general Newtonian gravity 
treatment in HJA, the integrals we obtain allow us to 
verify that a maximal uniform strain continues to yield 
the maximum quadrupole deformation in the Newtonian 
and relativistic no-Cowling cases. (UCB showed this to 
be true for an arbitrary equation of state in the Newto- 
nian Cowling approximation case; we are able to verify 
that it is true in the more general cases for each back- 
ground stellar model we consider.) 

We then apply our calculation to the standard case of 
quadrupoles supported by shearing the lattice of nuclei 
in the crust, as well as the cases where the quadrupole 
is supported by the hadron-quark mixed phase lattice 
in the core, or a crystalline color superconducting phase 
throughout a solid strange quark star. For the crustal 
quadrupoles, we calculate the shear modulus following 
HJA, using the equation of state (EOS) and composi- 
tion results of Douchin and Haensel [23| and the effective 
shear modulus calculated by Ogata and Ichimaru [24| . 
(There are recent improvements to the Ogata and Ichi- 
maru result [25l - f27| , but these reduce their shear modulus 
only by < 10%.) For the hadron-quark mixed phase, 
we use our recent calculations of the EOS and shear 
modulus for a variety of parameters. (We also con- 
sider the range of surface tensions for which the mixed 
phase is favored.) For crystalline quark matter, we use 
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the shear modulus calculated by Mannarelli, Rajagopal, 
and Sharma (29[, and the EOS given by Kurkela, Ro- 
matschke, and Vuorinen [30l |. 

In all cases, we use a breaking strain of 0.1, compara- 
ble to that calculated by Horowitz and Kadau [3l[ using 
molecular dynamics simulations. (Hoffman and Heyl |32j 
have recently obtained very similar values over more of 
parameter space.) This result is directly applicable to the 
crustal lattice (though see Chugunov and Horowitz [33[ 
for caveats). We also feel justified in applying it to the 
mixed phase and crystalline quark matter, since the pri- 
mary source of the high breaking strain appears to be 
the system's large pressure. But one can apply our re- 
sults to any breaking strain using the linear scaling of the 
maximum quadrupolc with breaking strain. 

In our general relativistic calculation, we use the rel- 
ativistic theory of elasticity given by Carter and Quin- 
tana (34j and placed in a more modern guise by Karlovini 
and Samuelsson (35| . However, all we need from it is 
the relativistic form of the clastic stress-energy tensor, 
which can be obtained by simple covariance arguments, 
as noted by Schumaker and Thorne 1301 . We also use the 
standard Thorne and Campolattaro [37| Regge- Wheeler 
gauge [38| formalism for perturbations of static relativis- 
tic stars, following Hinderer's recent calculation [3!| of 
the quadrupolc moment of a tidally deformed relativistic 
star (first discussed in Flanagan and Hinderer J!(|), and 
the classic calculation by Ipser [41 1. 

Even though we are interested in the gravitational ra- 
diation emitted by rotating stars, it is sufficient for us 
to calculate the static quadrupolc deformation. As dis- 
cussed by Ipser fi~0 | , and then proved for more general sit- 
uations by Thorne [42| . this static quadrupolc (obtained 
from the asymptotic form of the metric) can be inserted 
into the quadrupole formula to obtain the emitted grav- 
itational radiation in the fully relativistic, slow-motion 
limit. [This approximation has uncontrolled remainders 
of order (w/wr-) 2 , where w and uik are the star's angular 
velocity and its maximum — i.e., Kepler — angular veloc- 
ity, respectively. This ratio is < 10 -2 for the pulsars for 
which LIGO has been able to beat the spindown limit Q •] 

We shall generally show the gravitational constant G 
and speed of light c explicitly, though we shall take 
G = c = 1 in most of Sec. IIII1 only restoring them in our 
final expressions. The relativistic calculation was aided 
by use of the computer algebra system Maple and the as- 
sociated tensor manipulation package GRTensorII [43j . 
We used Mathematica 7 to perform numerical compu- 
tations. 

The paper is structured as follows: In Sec. HH we 
review UCB's formalism and extend it by introducing 
a Green function to compute the maximum Newtonian 
quadrupole deformation without making the Cowling ap- 
proximation. In Sec. IIII1 we further generalize to the 
fully relativistic case, and compare the various approx- 
imations for the maximum quadrupolc. In Sec. IIV1 
we show the maximum quadrupolcs for three different 
cases: first crustal quadrupoles, then hadron-quark hy- 



brid quadrupolcs, and finally solid strange quark star 
quadrupoles. We also describe the modifications to our 
formalism needed to treat solid strange quark stars. We 
discuss all these results in Sec. [Vj and summarize and 
conclude in Sec. IVII In the Appendix, we discuss uncer- 
tainties in the extent to which the mixed phase is favored 
in stable stars, especially when the total gravitational 
binding energy of the star is taken into account. 



II. NEWTONIAN CALCULATION OF THE 
MAXIMUM QUADRUPOLE 

We first demonstrate how to compute the maximum 
Newtonian quadrupole without making the Cowling ap- 
proximation. This provides a warm-up before we tackle 
the full relativistic case, and also allows us to verify some 
of the statements made by UCB and HJA. We use the 
basic formalism of UCB, modeling the star as nonrotat- 
ing, with the stress-energy tensor of a perfect fluid plus 
shear terms, and treating the shear contributions as a 
first-order perturbation of hydrostatic equilibrium. This 
pcrturbative treatment should be quite a good approx- 
imation: The maximum shear stress to energy density 
ratio we consider in the crustal and hybrid star cases is 
< 0.05% (and the maximum shear stress to pressure ra- 
tio is < 0.3%). And even in the case of solid strange 
quark stars, the maximum shear stress to energy den- 
sity ratio is still only at most ~ 0.2%. [We have already 
discussed the effects of rotation in the relativistic case, 
above; UCB note at the beginning of their Sec. 4 that 
rotation also only modifies the perturbative Newtonian 
results for the static deformations we and they consider 
at the 0([u/lo k ] 2 ) level] 

It is convenient to start by writing the quadrupole mo- 
ment in terms of the surface value of the perturbation to 
the star's Newtonian potential. We start from UCB's 
definition of 



5p(r)r A dr 



(1) 



(where the density perturbation Sp and all similar per- 
turbed quantities have only an I = m = 2 spherical har- 
monic component). [Note that this quadrupole moment 
differs by an overall constant from the one defined by 
Thorne Hi]— e.g., his Eq. (5.27a).] We then recall that 
the perturbed Poisson equation for the I = 2 part of the 
perturbed gravitational potential is 
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(i 



(A a <y#)(r) := ^[r 2 <5$'(r)]' - -j*$(r) = 4nG5p (2) 

(A2 is the I = 2 radial part of the Laplacian), with 
boundary conditions of 



5$(0) = 0, R6&(R) = -36§(R), 



(3) 



where R is the radial coordinate of the surface of the star. 
[See, e.g., Eqs. (2.15) and (2.16) in Q — their $22 is our 
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(5$. Note also that the primes denote derivatives with 
respect to r. Additionally, we shall continue to be incon- 
sistent with our inclusion of the functional dependence 
of quantities — e.g., Sp depends upon r, even though we 
do not always indicate this explicitly. We will eventually 
stop displaying <5$'s explicit functional dependence on r, 
for instance.] If we now substitute Eq. @ into Eq. (TTJ) 
and integrate by parts using the boundary conditions ([3]) , 
we obtain 



ST? 3 

Qr2 _ -——S$(R). 



(4) 



This sort of expression is more commonly seen in the 
rclativistic case, where it is necessary to obtain the 
quadrupolc in this manner by looking at the perturba- 
tion's asymptotic behavior — see the discussion in Sec. IIIII 
We now wish to obtain an equation for <5$ in terms of 
the shear stresses. We follow UCB in decomposing the 
perturbed stress tensor as [see their Eqs. (59) and (61)] 



Sr a b = SpYi m g ab + t rr Y lm (r a r b - e ab /2) + t r ^f ab 
+ t A (A a b+Yi m e ab /2). 



(5) 



Here 5p is the pressure perturbation, Y\ m is a spherical 
harmonic, f a is the radial unit vector, t rr , t r ±, and t\ are 
the components of the shear stresses, and g a b denotes the 
metric of flat, 3-dimensional Euclidean space. (Following 
UCB, we will generally write out I and m explicitly, even 
though we only consider I = m = 2 here.) Also [Eqs. (40) 
in UCB1, 



e-ab ■= 9ab - r a r b , 
f a b := 2rf( a V;,)Yi m //3, 

fl := y/l(l + l) = V6, 
Kb := r 2 V a V b Y lm /p 2 - 



fab/P- 



(6a) 
(6b) 

(6c) 
(6d) 



(We have corrected the dropped factor of (3 1 multiplying 
fab in UCB's definition of A a b — this was also noticed by 
HJA HI.) We also have 



We then project Eq. <j8j) along V b Yi m to express Sp in 
terms of the shear stresses t rr , t r ±, and along with p 
and (5$, giving 



Substituting this into Eq. (0), we thus obtain 



*A. (10) 



A 2 5$ — p 5$ 



9(r)' 



9(r) 



_f' /' f' f 

2 rr ^ r J_ n u rJ- 



1 1\ . 3 

p z 2 J T T 



(11) 



We now wish to obtain an integral expression for Q22 
that generalizes UCB's Eq. (64) to the case where we do 
not make the Cowling approximation. We shall do this 
by obtaining the Green function for the left-hand side of 
Eq. (fTTj) and then integrating by parts. We will be able 
to discard all of the boundary terms, since the stresses 
vanish at the star's surface (we assume that the shear 
modulus vanishes there) and the integrand vanishes at 
the star's center. We can obtain the Green function using 
the standard Sturm-Liouville expression in terms of the 
solutions of the homogeneous equation [e.g., Eq. (10.103) 
in Arfken and Weber [45| ]. We obtain the appropriate 
solution to the homogeneous equation numerically for a 
given background stellar model (EOS and mass). The 
equation for the Green function is [multiplying the left- 
hand side of Eq. (|lip by r 2 to improve its regularity] 



d_ 

dr 

= 8(r — f) 



or 



„ AirGr 2 . 
6+——p' 
9{r) 



Q(r, r) 
(12) 



tn 



^P^abi 



(7) 



where p is the shear modulus and o a b is the strain ten- 
sor. (This is a factor-of-2 correction to the expression in 
UCB, as noted in [H|.) Now, a convenient expression can 
be obtained from the perturbed equation of hydrostatic 
equilibrium 



V a <5r ab = 5pg{r)f b + pV b S^, 



(8) 



where V a denotes the flat-space covariant derivative, by 
substituting for Sp using the Poisson equation ^ and 
projecting along r b , yielding 



P 



4nG g(r) 



5& 



f b V a 5T ab 



9(r) 



1 

9(r) 



-op + t H — t rr t r j_ 

r r 



(9) 



with boundary conditions (at the star's center and sur- 
face) of 



G(0, f) = 0, Rd ± g(R, f) = -3Q(R, f), 



(13) 



where d\ denotes a partial derivative taken with respect 
to the first "slot" of the function. 

If we then write [using Eq. ((3]), the factor of r 2 from 
the Green function equation (|12[) , and the prefactor on 
the right-hand side of Eq. ([TT])] 



G N (r) :=-5R 3 r 2 g{R,r)/g{r), 



(14) 



we have 
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) 22 



G N (r) 



3 , 4 . r „ 



/3 



-G' N (r) - — Gjv(r) 



i rr + 



1 1\ , 3 /3 
p z 2 J r r 



1 1 

2 ~ Jp 



(15) 



G^(r)t A Ur. 



We have freely integrated by parts in obtaining the sec- 
ond expression, noting that the boundary terms are zero 
since Gjss(t) vanishes sufficiently rapidly as r — > and 
the stresses are zero at the surface of the star (since we 
assume that the shear modulus vanishes at the star's sur- 
face). 1 This reduces to UCB's Eq. (64) if we take the 
Cowling approximation 



G N {r) -+ r 4 /ff(r), 



(16) 



corresponding to dropping the second term on the left- 
hand side of Eq. (ITU . 

To obtain an analogue of the expression for the maxi- 
mum quadrupole given in Eq. (5) of Owen we note 
that UCB's argument about maximum uniform strain 
leading to the maximum quadrupole still holds here for 
the stars we consider, since the coefficients of the stress 
components in the integrand are all uniformly positive. 
(We have checked this numerically for each background 
stellar model we have considered.) The strain tensor 
components are 



G rr = (327r/15) 1/2 ff n 
a r± = (3/2)V2, 



a 



rr • 



(7a = 3oy 



(17a) 
(17b) 
(17c) 



in the case where the star is maximally (and uniformly) 
strained — see Eqs. (67) in UCB. The breaking strain <7 max 
is given by the von Mises expression, 



ab 



2ai 



(18) 



It thus corresponds to assuming that the lattice yields 
when it has stored a certain maximum energy density. 
We then have 



\Q 



max, N I 
22 



^ / p{r)[rG%(r)+iG' N (r)]dr. 



(19) 



1 We shall treat the case where the stresses do not vanish at the 
surface of the star when we consider solid strange quark stars in 
Sec. lIV Cl Also, note that HJA claim that UCB's expression does 
not include distributional contributions due to sudden changes 
in the shear modulus. This is not the case — these are included 
due to UCB's integration by parts (cf. the definition of the dis- 
tributional derivative). All that the UCB derivation requires is, 
e.g., that the shear modulus vanish outside of the crust, not that 
it do so continuously. 



This reduces to Eq. (5) in Owen [12| if we use the Cowling 
approximation [Eq. (|16[) ]. 

Note that there is no direct contribution from p' to G'^ 
in the no-Cowling case, despite what one might expect 
from Eq. (|T2"T) : Writing Q(r) := Q(R,r) for notational 
simplicity, the p' contribution from 

§"(r) = (2/r)£'(r) + [6/r 2 + ^Gp'{r)/g{r)]g{r) (20) 

is exactly canceled by one from 

g"(r) = 6Gm(r)/r 4 - ^Gp{r)/r + 4irGp'(r) (21) 



G%(r) = -5R 3 r 2 ig"(r)/g(r)-g(r)g"(r)/{g(r)¥ 
+ {terms with no p'}]. 



(22) 



However, there is a direct contribution from p' to G'^ (via 
g") if we make the Cowling approximation [Eq. ([TB")) ]. We 
shall see that this leads to a significant difference in the 
resulting contributions to the quadrupole moment from 
regions of the star surrounding a sudden change in den- 
sity (e.g., near the crust-core interface, which will be rel- 
evant for the quadrupolcs supported by crustal elasticity 
considered by UCB and others). 

Numerically, we compute 67 at using the standard ex- 
pression for the Green function in terms of the two in- 
dependent solutions to the homogeneous equation [see, 
e.g., Eq. (10.103) in Arfken and Weber [ii] ]. Since we 
are solely interested in the Green function evaluated at 
the star's surface, we can eliminate one of the homoge- 
neous solutions using the boundary conditions there, and 
only consider the homogeneous solution that is regular 
at the origin, which we call F . In terms of F, the Green 
function is given by 



g(R,r) = - 



F(r) 



3RF(R) + R 2 F'(R) ' 



(23) 



We thus solve CnF = [with the operator Cn given by 
Eq. (|T2"j) ] with the boundary conditions -F(ro) = 1 and 
F'(ro) = 2/ro, where rg is the small inner radius used 
in the solution of the OV equations, as discussed at the 
end of Sec. IIIII [These boundary conditions come from 
regularity at the origin, which implies that F(r) = 0(r 2 ) 
there.] 

Our Green function method for obtaining the max- 
imum quadrupole numerically may seem more compli- 
cated than existing methods because it introduces extra 
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steps. However this method is ideal for showing that 
maximum stress gives the maximum quadrupole and for 
seeing how much stresses at different radii contribute to 
the total quadrupole. It also appears to be the simplest 
way of dealing with any potential distributional contri- 
butions from the derivatives of the shear modulus, since 
they are automatically taken care of by the integration 
by parts. 



III. GENERAL RELATIVISTIC CALCULATION 
OF THE MAXIMUM QUADRUPOLE 

Here we compute the maximum quadrupole moment in 
general relativity, using the Regge- Wheeler gauge [HI rel- 
ativistic stellar perturbation theory developed by Thorne 
and Campolattaro [I?} , as in the similar calculation of the 
tidal Love number of a relativistic star by Hinderer (3t| . 
We start by writing down the line element corresponding 
to a static, even-parity, I = 2 first-order perturbation of 
a static, spherical, relativistic star in the Regge- Wheeler 
gauge [cf. Eq. (14) in Hinderer [U ]: 

ds 2 = -[1 + H (r)Y lm ]f(r)dt 2 + [1 + H 2 (r)Y lm ]h(r)dr 2 



[l + K(r)Y lm ]r 2 (d6 2 



(24) 



Here we have used the notation of Wald [46[ for the back- 
ground, so that / and h are the standard Schwarzschild 
functions for the unperturbed star, with / = e 2 *, where 



<j>'[r) 



m(r) + 4irr 3 p 



[r — 2m(r)] 
with <j>{R) = log(l - 2M/R)/2, and 

2m(r) 



h{r) 



In these expressions, 



m(r) := 47r / p(r)r 2 dr 



(25) 



(26) 



(27) 



Also, recall that we write our spherical harmonics in 
terms of I and m, following UCB, even though we special- 
ized to I = m = 2, and that we are now taking G = c = 1. 

The metric perturbation is determined by Hq, H 2 , and 
K, which here are sourced by the perturbation to the 
star's stress-energy tensor. The appropriate stress-energy 
tensor can be obtained directly from the standard New- 
tonian expression ([5]) by simple covariance arguments, 
as in Schumaker and Thorne [36| . or from the detailed 
relativistic elasticity theory of Carter and Quintana [IH 
[see their Eq. (6.19); this is also given in Eq. (128) of 
Karlovini and Samuelsson [35[ ]. All we really need for 
our purposes is to note that the shear contribution is 
tracefree with respect to the background metric, so that 



we can use the obvious covariant generalization of the 
decomposition given by UCB, 2 yielding 

ST ab = [SptJb + Sp(g a b + iJb) - t rr (f a f b - q a b/2)]Yi m 
- t r ±f ab - t A (A ab + h 1/2 Y lm q ab /2), 

(28) 

with the full stress-energy tensor given by 

T a b = ptj> + p(S a b + U b ) + ST a b . (29) 

Here, indices now run over all four spacetime dimen- 
sions and g a b denotes the background (spacetime) met- 
ric (which we use to raise and lower indices). Addition- 
ally, we have introduced the background temporal and 
radial unit vectors t a and r a ; q a b is the induced met- 
ric on the unit 2-sphere; /„& := 2rf( a \?b)Yi m //3] and 

A Qb := r 2 /i 1 /2 Va V b l, m //? 2 + fab/ p. Here f„ and V Q 
now have their curved-space meanings. 

Our A a b differs from the Newtonian A a b [from UCB, 
given in our Eq. ([6d| ] due to the insertion of h 1 / 2 . This 
insertion is necessary for A ab to be transverse and orthog- 
onal to fab (with respect to the background spacetime 
metric). The same logic leads to the introduction of the 
factor of h 1 / 2 multiplying q a b in the t\ term in Eq. (|28[k 
it is there so that the t\ term is orthogonal to the t rr 
term. We have used UCB's convention for the relative 
sign between the perfect fluid and shear portions of the 
stress-energy tensor, though we have reversed the overall 
sign. (However, we used the UCB convention proper in 
Sec. [TTJ) The factor of /i 1 / 2 in the coefficient of t\ leads 
to a factor of hr 1 in the strain a\ that corresponds to 
the von Mises breaking strain (TIB)) . We thus replace the 
Newtonian Eq. (|17cp with 



(TA = 3a rr /h, 



(30) 



leaving Eqs. Q17ap and (|17b[) unchanged. 

One can now obtain an equation for Hp from the per- 
turbed Einstein equations, as in Ipser 41(. (The other 
two metric functions, H 2 and K, can be expressed in 
terms of H ; these expressions are given by Ipser.) The 
concordance for notation is v = 2<fi, e v = /, A = 2ip, 
e x =h,pi = -Sp, pi = -5p, <£ 2 = t rr , Qi = h 1/2 W/?, 
and & = h 1 ' 2 t\/P' 2 . Additionally, Ipser's H is the 
negative of ours. The relevant result is given in Ipser's 



2 Of course, this assumes that it is possible to obtain any sym- 
metric tracefree tensor from the detailed relativistic expression, 
but — as would be expected (and can easily be seen from the 
expressions) — this is indeed the case, at least if one only works 
to first order in the perturbation, as we do here. Also, it is 
instructive to note that we do not need to know the specifics 
of the matter displacements that generate the quadrupoles we 
consider, only that there is a tracefree contribution to the star's 
stress-energy tensor whose maximum value is given by the ma- 
terial's shear modulus and von Mises breaking strain. 
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Eqs. (27)-(28), and is 

Hg + + <f>' - H + V{r)H = 8TTh^ 2 S(r), (31) 
where 
V(r) 



20" + 24>' ( - - $ - A + ^ - (32) 
1 r r r £ 



and 
Sir) :-- 



- h l /\8p + Sp t rr ) + 2 j(3 - r0')^ + 
+ [rV + r«^(5 - r<^) + rtj}' - f3 2 h/2 + l]^r 



/3 2 



2 J.'^A 

h 1 / 2 (Sp + 5p)+S [t] (r) 7 



(33) 

in our notation. Here we have defined ip := (1/2) log h 
and written for the contributions from shear stresses. 
(The "=:" notation implies that the quantity being de- 
fined is on the right-hand side of the equality.) 

We now wish to eliminate Sp and Sp in favor of the 
shear stresses, as in the Newtonian calculation. We use 
the same projections of stress-energy conservation as in 
the Newtonian case (projecting onto the quantities de- 
fined by the background spacetime, for simplicity) along 
with the Oppenheimer-Volkov (OV) equations, giving 



and 

s p = —^-(p+p) - 



Ph 1 / 2 



(34) 

[(3 + r<t>')t rl _ + < J 
(35) 

Using the second expression to substitute for Sp' in the 
first, we have 



S P + Sp = U^(p'+p') 



f3 2 h 



t + -t' 



— + </>' + r<p" - V'[3 + r<f>'] t r± 
r 



+ (4 + rfe'-^]K± + <± 
x h 1/2 (^'t A + t' A 



1 



(36) 



The equation for Ho thus becomes 

'2 



(£ GR H )(r) := H 



<t>' - ^) #o 



V(r) - 4irh- 



H 



8nh 1 ' a [h 1 ' 2 S [S p,Sp]( r )/<t>'+ S lt]( r )]- 



(37) 



[P(r) and <Sm(r) are given in Eqs. (|32|) and (|33|) . respec- 
tively] As expected, this reduces to Eq. (fTT|) in the New- 
tonian limit [where we have Ho — > 2(5$ and <f>' — > g(r)]. 

We now want to write the equation for H in Sturm- 
Liouvillc form in order to obtain its Green function easily 
To do this, we note that the appropriate "integrating 
factor" (for the first two terms) is r 2 (f /h) 1 / 2 , which gives 



lr 2 (f/hy/ 2 H>}> + r 2 (f/h)V 2 



V(r) - 4tt/i 



P'+P' 



Ho 



= 87rr 2 / 1/2 [^ 1/2 5 [5p , 5p] (r)/0' + S [t] (r)}. (38) 

We also need the boundary conditions, which are given 
by matching Ho onto a vacuum solution at the surface of 
the star. The vacuum solution that is regular at infinity 
is given by Eq. (20) in Hinderer [39j with c 2 = 0, viz., 



Ho(R) = ci 



I- 1 



C 2 /2 + 3C-3 
l~~C 



where we have evaluated this at the star's surface (r 
and defined the star's compactness 

C := 2GM/Rc 2 



(39) 



R) 



(40) 



(now returning to showing factors of G and c explic- 
itly). We require that Ho and H' be continuous at the 
star's surface. The value of c\ obtained from this match- 
ing of the internal and external solutions gives us the 
quadrupole moment. If we use the quadrupolc moment 
amplitude that reduces to the UCB integral [given in our 
Eq. ([T])] in the Newtonian limit, we have 



Q-i 



G 2 M 3 Cl 



(41) 



[This expression comes from inserting a pure I = m = 2 
density perturbation into Eq. (2) in Hinderer [39| . con- 
tracting the free indices with unit position vectors, per- 
forming the angular integral, for which the expressions 
in Thorne [42[ are useful, and noting that the result is 
(87r/15)l22 times our Eq. (TTJ. The given result then fol- 
lows immediately from Hinderer's Eqs. (7), (9), and (22); 
we reverse the overall sign since we have reversed the 
UCB sign convention for the stress-energy tensor.] 
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We then have a Green function for Q22 of 



where 



S GR (R,r) 



2GMV 
Rc 2 ) 
U(r) 



(42) 



c 2 R 2 [U'(R)H a (R) -U{R)H' (R)} 



(including the overall factor of 8nG/c 4 that multiplies 
the source). Here U is given by C GR U = [£ G r 1S 
given by Eq. (|37[) ]. with boundary conditions U{tq) = 1 
and U'{r ) = 2/r . [Compare Eq. (10.103) in Arfken 
and Weber (4f| , as well as our Newtonian version above.] 
Here Hq(R) and H' (R) are given by the boundary con- 
ditions (|39p with ci — > 1. [One obtains this expression by 
first computing the Green function for Hq(R) following 
Arfken and Weber, then dividing through by the quan- 
tity in brackets in Eq. Q39[) to obtain ci, and finally us- 
ing Eq. ([4"Tjl to obtain Q22- We have also noted that 
Iff -> h -> 1/(1 - 2G1\I/Rc 2 ) at the star's surface.] 
We thus define, for notational simplicity, two rclativistic 
generalizations of GAr(r): One, 



G GR (r) := 



r 2 {fhf/ 2 g GK (R,r) 

4>> 



(43) 



for the contributions from S\s P .Sp] 1 an d one j 

G GR (r) :=r 2 f^ 2 g GR (R,r), (44) 

for the contributions from <Sr t i . 

With these definitions, the integral expression for the 
quadrupolc in terms of the stresses and the structure of 
the background star is 



22 = / [G GK {r)S [SpiSp ] (r) + G GR {r)S [t] (r)] dr 



(Crrtrr + C t ±t r ± + C\t\)dr, 



(45) 



C rr := + 4>^j G GR (r) - |G' GK (r) - h^ 2 G GR (r), 

(46a) 



C r _ 



/3k 1 / 2 



2 + r(£+£) 
G grM + ^T72 G grW 



/3/i 



TT^grM + l_M<5 GR (r) - 2^G GR (r), 



(46b) 



-G GR (r) (46c) 



2r^' (3 - r<£') + 2rV' - P 2 h + 2 
+ ^ 

Jg2 GRv J> 

and we have integrated by parts twice to obtain the sec- 
ond equality in Eq. (|45p. using the same argument as in 
our Newtonian calculation. 



We now look at the maximum quadrupole. This is 
still given by the uniformly maximally strained case: We 
have checked numerically that the coefficients of the three 
stress terms are always negative for all the background 
stars we consider. We thus have a maximum quadrupolc 
given by inserting Eqs. ©, (|T7al) . (|17b|) . and ([30) into 
Eq. (|S5)1. yielding 



i^max.GRi 
I V22 



where 



327T f R 

Wo " {r) 



-(h 1 / 2 -l)-2<p> 



G GR (r) 



) + [ 3 - + G' GR (r) + ^G GR (r) + Q strcss jdr, 



(47) 



s 



rW3)-^-l + 1/2 + 1 



(48) 



is the contribution from the stresses' own gravity. We 
have split it off both for ease of notation and because it 
is negligible except for the most massive compact stars, 
as illustrated below. The contributions from the density 
and pressure perturbations are so much larger due to the 
factor of 1/0' present in G G r [cf. Eqs. (|43]) and flU}]. It 
is easy to see that Eq. (|4"T) reduces to Eq. (fH?)) in the 



Newtonian limit, where h 1, and we can neglect the 
contributions involving cf>' , 7//, and Q stross . 

We now show how the different maximal-strain Q22 
Green functions [given by the integrands in Eqs. f| 19[) 
and (|4"T|) without the factors of /j] vary with EOS, as well 
as with the mass of the star for a given EOS. This gives an 
indication of how much difference the various approxima- 
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r(km) r(km) 



FIG. 1: The Q22 integrands (without the factor of /icr max ) for 
the SLy EOS and an 0.500 Mq star with a compactness of 
0.12. 



FIG. 3: The Q22 integrands (without the factor of fia ma , x ) 
for the SLy EOS and a maximum mass, 2.05 Mq star, with a 
compactness of 0.60. 




r(km) 



FIG. 2: The Q22 integrands (without the factor of /JCTmax) for 
the SLy EOS and a 1.40 Mq star with a compactness of 0.35. 



tions make in different situations. We start with the uni- 
fied SLy EOS (HI, obtained by Haensel and Potekhinjp] 
(using the table provided by the IofFc group [H| at [49j|). 
which is a standard choice for making predictions about 
crustal quadrupoles (e.g., in HJA and our Sec. IIV A[) . 
Here we illustrate the changes in the Green functions 
with mass for stars with masses ranging from 0.5 Mq to 
the EOS's maximum mass of 2.05 Mq; see Figs. [TJ [3J and 
[3J (All three Green functions agree extremely closely 
for stars around the EOS's minimum mass of 0.094 Mq, 
so we do not show this case, particularly because such 
low-mass neutron stars are of unclear astrophysical rele- 
vance.) These stars' compactnesses [defined in Eq. (|4"0J)] 
range from 0.12 to 0.6. Note that Fig. [3] has a different 
vertical scale than the other two plots, due to the sup- 
pression of the quadrupole for massive, compact stars 
(discussed below). 

We illustrate the ratios of the various Q22 Green func- 
tions to the Newtonian Cowling approximation one for 



FIG. 4: Ratios of the Q22 integrands with the Newtonian 
Cowling approximation integrand for the Hyl EOS and a 
maximum mass, 2.06 Mq star, with a compactness of 0.49. 
Note that the top and bottom plots have completely separate 
vertical axis scalings. 



the maximum mass (2.05 Mq) hybrid star using the Hyl 
EOS (see Table I in [28| ) in Fig. g] 3 We see the overesti- 
mate of the Newtonian no Cowling approximation calcu- 
lation for perturbations in the core, particularly com- 
pared with the general relativistic (GR) version, and 
also see the overestimate of the Newtonian Cowling ap- 



As discussed in 1281 . for our low-density EOS, we use the same 
combination of the Baym, Pethick, and Sutherland (BPS) 5y| 
EOS for n B < 0.001 fm~ 3 and the Negele and Vautherin |51[ 
EOS for 0.001 fm~ 3 < n B < 0.08 fm" 3 used by Lattimer and 
Prakash [52l | . These were obtained from the table provided by 
Kurkela et al. [53| at |54| . Bulk quantities of hybrid stars such 
as the mass and quadrupole moment (from core deformations) 
do not depend much on the precise choice of low-density EOS. 
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proximation version for crustal perturbations. (We do 
not make some sort of similar plot for the solid strange 
quark star case, since the expressions for the maximum 
quadrupole in this case end up being rather different than 
the intcgrated-by-parts ones presented in the previous 
sections, as we shall see in Sec. lIVCl ) 

In all these cases, we compute the stellar background 
fully rclativistically, using the OV equations and iden- 
tifying the OV equations' Schwarzschild radial coordi- 
nate with the Newtonian radial coordinate when neces- 
sary. We have used the enthalpy form of the OV equa- 
tions given by Lindblom [55| and implemented the inner 
boundary condition by taking the star to have an in- 
ner core of radius tq = 100 cm, whose mass is given by 
(4/3)7rrQ£o, where eo is the energy density correspond- 
ing to the central enthalpy that parametrizes the solu- 
tion. (The spike near the origin seen in the bottom 
plot in Fig. [4] is due to this implementation of the in- 
ner boundary condition and has a negligible effect on the 
computed maximum quadrupoles.) In all cases, we have 
used Mathematica 7's default methods to solve the dif- 
ferential equations, find roots, etc. We have computed as 
many derivatives as possible analytically, to aid numeri- 
cal accuracy, e.g., using the OV equations to substitute 
for derivatives of the pressure, and also using the Green 
function equations to express second derivatives of the 
Green functions in terms of the functions themselves and 
their first derivatives. 



IV. RESULTS 
A. Maximum Q22 for crustal deformations 

Here we consider the maximum quadrupoles from clas- 
tic deformations of a nonaccreted crust in three possi- 
ble situations, following HJA. In particular, we use the 
SLy EOS (as do HJA, though they do not refer to it 
by that name) and impose two comparison crustal thick- 
nesses to ascertain how much this affects the maximum 
quadrupole. Here we use the same rough model for the 
crust's shear modulus used by HJA. Wc also consider 
the more detailed model for the shear modulus obtained 
using the crustal composition provided by Douchin and 
Haensel (23[. Here the crust's thickness is fixed to the 
value given in that work. In this case, we also consider a 
different high-density EOS that yields much less compact 
stars with larger crusts. 

Specifically, the two comparison crustal thicknesses are 
given by taking the base of the crust to occur at densi- 
ties of 2.1 x 10 14 g cur 3 (thick crust, for comparison with 
UCB) or 1.6x 10 14 _g cm -3 (thin crust, following a sugges- 
tion by Haensel |5g|), while Douchin and Haensel place 
the bottom of the crust at a density of 1.28x 10 14 g cm -3 . 
For the two comparison cases, we take the shear modulus 
to be 10 16 cm 2 s -2 times the star's density (in g cm -3 ). 
As illustrated in HJA's Fig. 2, this is an underestimate 
of < 50%, except at the very extremes of the density 



range considered. 4 We plot the quadrupole moment 
and ellipticity for these two cases for masses between 
~ 1.2 Mq (around the minimum observed neutron star 
mass — see [571 ]) and the SLy EOS's maximum mass of 
2.05 M Q in Fig. 



1 


1 1 1 1 1 1 1 


| 1 | 1 | 1 | 1 | 1 






















■ ■ ■ Newt. Cowling (thick) 






Newt, no Cowling (thick) 






GR (total, thick) 






— Newt. Cowling (thin) 






Newt, no Cowling (thin) 






■■- GR (total, thin) 




■ 1 


, 1,1,1, 


1,1,1,1,1,- 



10 



M (solar masses) 

FIG. 5: The Newtonian Cowling, Newtonian no Cowling, and 
full relativistic (including stress contributions) values for the 
maximum quadrupole deformations (and fiducial ellipticity) 
due to crustal stresses versus mass for two choices of crustal 
thickness. These are computed using the SLy EOS with the 
rough HJA recipe for the shear modulus and a breaking strain 
of 0.1. 

In addition to the quadrup ole mo ments, we also show 
the fiducial ellipticity egd = ^/87r/15Q22/ I zz [e.g., Eq. (2) 
of [IH ]• Here I zz is the star's principal moment of inertia, 
for which we use the fiducial value of I zz = 10 38 kg m 2 = 
10 45 g cm 2 used in the LIGO/Virgo papers rather than 
the true value for a given mass and EOS, which can be 
greater by a factor of a few. We do this for easy compar- 
ison with the observational papers, since they frequently 
quote results in terms of this fiducial ellipticity instead 
of the quadrupole moment, which is the quantity truly 
measured. 

Nota bene (N.B.): We present these fiducial ellipticities 
only for comparison with LIGO/Virgo results, not to give 
any indication of the size of the deformation. While the 
true ellipticity gives a measure of the size of the deforma- 
tion in the Newtonian case (up to ambiguities due to the 
fact that the true density distribution is nonuniform), it 
does not do so in any obvious way in the relativistic case. 
If one wished to know the true, relativistic size of the de- 
formation (or the shape of the deformed star's surface), 
one would need to calculate it using the same perturba- 
tive setup used to obtain the quadrupole moment. (This 
would also necessitate using a detailed relativistic theory 



4 Note that Fig. 3 in HJA is not in agreement with their Fig. 2. 
When we reproduce those figures, we find that the ratio fi/p is 
considerably closer to 10 16 cm 2 s — 2 over all the density range 
than the trace shown in HJA's Fig. 3, so their approximation is 
better than it would appear from that figure. 
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of elasticity to relate the stresses to the matter displace- 
ments.) 

In the more detailed case, we use the HJA version 
of the Ogata and Ichimaru [HI shear modulus, com- 
bined with the Douchin and Haensel [23] results for the 
crust's composition. This is [correcting a typo in HJA's 
Eq. (20)], 



Mcff 



0.1194 



4-7T 



1/3 



A 



4/3 

n b ) (Ze) 2 , (49) 



where X n is the fraction of neutrons outside of nuclei, A 
and Z are the atomic and proton number of the nuclei, 
respectively, n& is the baryon number density, and e is 
the fundamental charge. 

Since HJA's study, there have been a few improve- 
ments to the Ogata and Ichimaru result: Horowitz and 
Hughto [25[ have computed the effects of charge screen- 
ing, finding a ~ 7% reduction in the shear modulus. 
Baiko [26[ has also considered a relativistic model for 
the electron polarizability and arrived at similar conclu- 
sions. Indeed, Baiko's results suggest that screening will 
yield an even smaller correction in the innermost portion 
of the crust, where the shear modulus is the largest, and 
the electrons are the most relativistic, with a relativity 
parameter over an order of magnitude larger than the 
largest Baiko considers. (However, the ion charge num- 
bers are also almost always somewhat greater than the 
largest Baiko considers, particularly at the very inner- 
most portion of the crust, which will tend to increase the 
effect.) 

Baiko [13 has also recently computed quantum cor- 
rections, and finds that they reduce the shear modulus 
by up to ~ 18% in some regimes. However, in our case, 
the reduction will be much smaller, based on the scal- 
ing of p 1//6 / \Z A 2 / 3 ) given near the end of Baiko's Sec. 6. 
Even though our densities are over an order of magnitude 
greater, the nuclei we consider are also over an order of 
magnitude more massive than the 12 C composition Baiko 
considers, so the quantum mechanical effects end up be- 
ing reduced by about an order of magnitude from the 
number Baiko quotes. We thus use the same Ogata and 
Ichimaru result used by HJA, noting that the resulting 
quadrupoles might be reduced by less than 10% due to 
charge screening and quantum effects — an error which 
is small compared to other uncertainties, such as crust 
thickness and the composition of dense matter. Indeed, 
there is a factor of ~ 2 uncertainty in the shear modulus 
due to angle averaging (even disregarding whether the 
implicit assumption of a polycrystalline structure for the 
crust is warranted): As shown by Hill (5^|, the Voigt av- 
erage used by Ogata and Ichimaru is an upper bound on 
the true shear modulus of a polycrystal. A lower bound 
is given by the Reuss average (also discussed in Hill [58| ) , 
for which the prefactor in Eq. would be 0.05106. 

Note that there would be even further corrections to 
the shear modulus due to pasta phases (see [H|), but 
such phases are not present in the Douchin and Haensel 
model |23| . We also note that the Douchin and Haensel 



results only include the very innermost portion of the 
outer crust. However, this lack of coverage has a negli- 
gible effect on the- final results for the quadrupoles, since- 
the neglected region has at most half the radial extent 
of the inner crust and the shear modulus in this region 
is orders of magnitude below its maximum values at the 
bottom of the inner crust. We have checked this explicitly 
using the detailed calculations of the outer crust compo- 
sition due to Riister, Hempel, and Schaffncr-Biclich [6C 
available at 



61|. 



We plot the maximum quadrupole and ellipticity in 
the three approximations for the detailed shear modu- 
lus model in Fig. [51 Here we show these for the SLy 
EOS proper, and also for a high-density EOS that yields 
much less compact stars (and a crust that is ~ 2 times as 
thick), and thus larger maximum quadrupoles. For the 
latter EOS, we have chosen (for simplicity) the LKR1 
hybrid EOS from [28j — the maximum compactnesses for 
the two EOSs are 0.60 (SLy) and 0.43 (LKR1). (We shall 
show the much larger quadrupoles that could be sup- 
ported by the mixed phase in the core for the LKR1 EOS 
in Sec. IIVB1 but here just show the crustal quadrupoles 
using the Douchin and Haensel model for the crust.) 



Newt. Cowling (LKR1, crust) 
-■ Newt, no Cowling (LKR1, crust) 
- — GR (total, LKR1, crust) 
■ ■ ■ Newt. Cowling (SLy) 

Newt, no Cowling (SLy) 

GR (total, SLy) 



1.5 1.6 1.7 

M (solar masses) 



1.9 2 2.1 



FIG. 6: The Newtonian Cowling, Newtonian no Cowling, and 
full relativistic (including stress contributions) values for the 
maximum quadrupole deformations (and fiducial ellipticity) 
due to crustal stresses versus mass, for the SLy EOS with the 
detailed Douchin and Haensel + Ogata and Ichimaru model 
for the shear modulus and a breaking strain of 0.1, plus the 
crustal quadrupoles for the LKR1 EOS with the same crustal 
model. 

In all of these crustal results, in addition to the ex- 
pected relativistic suppression of the quadrupole (which 
becomes quite dramatic for compact, high-mass stars), 
we also find that the Newtonian Cowling approximation 
slightly overestimates the quadrupole (by ~ 25-50%), as 
observed by HJA (though they found the overestimate to 
be considerably greater, around a factor of at least a few). 
This overestimate is due to the cancellation of contribu- 
tions from p' when one drops the Cowling approximation 
(see the discussion at the end of Sec. [H]). The overall de- 
crease in the maximum crustal quadrupole with mass is 
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due primarily to the fact that the crust thins by a fac- 
tor of — 4 (SLy) or ~ 2 (LKR1) in going from a 1 M Q 
star to the maximum mass star, though the quadrupole 
itself receives even further suppressions with mass due to 
relativistic effects and an increased gravitational field. 



B. Maximum Q22 for hybrid stars 




M (solar masses) 



FIG. 7: The Newtonian Cowling, Newtonian no Cowling, and 
full relativistic (including stress contributions) values for the 
maximum quadrupole deformations (and fiducial ellipticity) 
of hybrid stars versus mass, using the Hyl EOS with a surface 
tension of a = 80 MeV fm~ and a breaking strain of 0.1. 

Here we display the maximum quadrupole deforma- 
tions as a function of stellar mass for each of the hy- 
brid EOS parameter sets considered in [HI . We start by 
showing these values calculated in the various approx- 
imations using the Hyl EOS (with a surface tension of 
g = 80 MeV fm~ 2 ; see Table I in [H) in Fig. El and then 
restrict our attention to the relativistic results. (The re- 
lation between the results of the different approximations 
is roughly the same for all the hybrid EOSs we consider.) 
Here the maximum quadrupoles increase with mass, since 
the volume of mixed phase increases with mass, and this 
is more than enough to offset the suppressions due to 
relativity and the increased gravitational field. 

We also show how the maximum relativistic 
quadrupole varies with the surface tension for the Hyl 
EOS in Fig. [5J The slightly larger quadrupoles for lower 
surface tensions at low masses are expected, due to a 
slightly larger shear modulus at low pressures for lower 
surface tensions — see Fig. 10 in [28j]. In fact, despite the 
greater than an order of magnitude differences in high- 
pressure shear modulus for the Hyl EOS in going from a 
surface tension of 20 MeV fm~ 2 to one of 80 MeV fnV 2 
(see Fig. 10 in [IH), the differences in the resulting maxi- 
mum quadrupoles are at most a factor of a few (for large 
masses). This is not unexpected: These quantities are 
dominated by the portions of the mixed phase further 
out in the star, where the shear moduli have a much 
weaker dependence on the surface tension. (Addition- 
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FIG. 8: The full relativistic maximum quadrupole deforma- 
tions (and fiducial ellipticity) of hybrid stars versus mass, 
using the Hyl EOS with various surface tensions, <r, and a 
breaking strain of 0.1. 



ally, the fact that larger surface tensions lead to smaller 
shear moduli at low pressures helps to minimize the ef- 
fect, though the maximum quadrupoles still increase with 
increasing surface tension for high masses, as expected.) 
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FIG. 9: The full relativistic maximum quadrupole deforma- 
tions (and fiducial ellipticity) of hybrid stars versus mass, 
using various EOSs from Table I in [2||, all with a surface 
tension of a — 80 MeV fm -2 and a breaking strain of 0.1. 

Finally, we show the maximum quadrupoles for dif- 
ferent hybrid EOSs in Fig. |H1 (Note that these curves 
start somewhat above the minimum masses for which 
the mixed phase is present, since we are mostly inter- 
ested in the significantly larger maximum quadrupoles 
possible for larger masses.) The considerable differences 
are due primarily to the substantial variations in the ex- 
tent of the mixed phase in stable stars with EOS param- 
eters as well as the EOS dependence of the stars' com- 
pactnesses (see Table I in [28||), not to variations in the 
magnitude of the shear modulus for a g iven quark mat- 
ter fraction (compared in Fig. 12 in [28J). In particular, 
the LKRl EOS produces stars with a very large region 
of mixed phase — up to 72.6% of the star's radius — and 
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a (relatively) small maximum compactness — only 0.433. 
(Note that our quadrupole curve for the LKR1 EOS ends 
slightly short of the EOS's maximum mass of 1.955 Mq, 
only going to 1.948 Mq, due to problems with the nu- 
merics.) 
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FIG. 10: The full relativistic maximum quadrupole deforma- 
tions (and fiducial ellipticity) of hybrid stars versus mass, for 
the LKR1 and Hyl' EOSs from [28| and a variety of surface 
tensions. Here, "plain" denotes no additions, while "full" de- 
notes both blob and lattice additions. For the full cases, the 
range of masses for which the mixed phase is favored by the 
binding energy argument (see the Appendix) is denoted by 
thick solid lines. 

N.B.: These maximum quadrupoles may all be overly 
optimistic. First, as was discussed in Sec. IIV A[ the av- 
eraging used to obtain the effective shear modulus only 
gives an upper bound on the true shear modulus of a 
polycrystal. (We do not quote results for the Rcuss 
lower bound here, since it is only straightforward to ob- 
tain for the three-dimensional droplet phases. However, 
we shall note that preliminary investigations, using the 
Reuss bound for the droplet phases, and the Voigt bound 
for the rest, give reductions in the maximum quadrupoles 
of up to ~ 5 for lower masses.) 

Second, as mentioned in Sec. II C of [28|], for large 
surface tensions the mixed phase is disfavored by local 
energy arguments, compared to the pure hadronic phase. 
It could still be possible for the mixed phase to be favored 
by global energy arguments. However, we find that even 
considering global energy arguments, the mixed phase is 
only favored for higher-mass stars with low surface ten- 
sions and the two most extreme EOSs from [28| (consid- 
ering surface tensions high enough that the entire lattice 
is stable), viz., LKR1 and Hyl'. (Though note that re- 
cent calculations place the surface tension on the low 
side — see [62| for the latest results.) But since this con- 
clusion depends rather sensitively on the energy density 
and pressure of the lattice, and our calculation in [28( 
does not include some contributions that would lower the 
lattice's energy, it is likely on the pessimistic side. It thus 
seems reasonable to take the results presented in Fig. [5] 
as upper bounds on the maximum quadrupole (and the 



range of masses for which the mixed phase is favored), 
and the binding energy argument results as rough lower 
bounds. In Fig. [TU1 we show the maximum quadrupoles 
for smaller surface tensions for those two EOSs, and in- 
dicate the masses for which the mixed phase is favored 
by the binding energy argument. We discuss this issue 
in more detail in the Appendix. 

Note that LIGO's current upper limits on fiducial el- 
lipticity in the most interestingcases (the Crab pulsar, 
PSR J0537-6910, and Cas A) [1,0] are - 1CT 4 , corre- 
sponding to a quadrupole moment of ~ 10 41 g cm 2 . The 
first hybrid star estimate by Owen fl2l ] was an order of 
magnitude lower. Thus our new results here show that 
current LIGO upper limits are interesting not only for 
quark stars but also for hybrid stars, at least high-mass 
ones. Indeed, the most extreme case we consider, the 
LKR1 EOS with high surface tensions, gives maximum 
quadrupoles of a few x 10 42 g cm 2 , which are above and 
therefore relevant to the limits set by Virgo for the Vela 
pulsar @- 



C. Maximum Q22 for crystalline color 
superconducting quark stars 

Here we consider stars made of crystalline color su- 
perconducting quark matter, for which the shear mod- 
ulus has been estimated by Mannarelli, Rajagopal, and 
Sharma [Hf. 5 [See Eq. (1) in Haskell et al. [H for the ex- 
pression in cgs units.] Such stars have also been treated 
(with varying degrees of sophistication) by Haskell et 
al. [13], Lin [3], and Knippel and Sedrakian (l9| . How- 
ever, only Lin considers the case of a solid quark star, as 
we will do here, and does so using quite a rough model. 
(The others consider crystalline color superconducting 
cores in hybrid stars.) 

Since strange quark stars have a nonzero surface 
density — and solid quark stars have a nonzero surface 
shear modulus, with the standard density-independent 
treatment of the superconducting gap — we have to make 
some changes to our previously obtained expressions in 
order to treat them. 

First, the outer boundary condition changes. The po- 
tential (in the Newtonian case) and metric perturbation 
(in the GR case) are no longer continuous at the star's 
surface, due to the presence of p' in both equations [see 
Eqs. ([T^]) and ([55]) ]. As discussed in Hinderer et al. [631 ] 
(following Damour and Nagar [64j]). one can obtain the 
distributional contribution to the boundary conditions 
using the usual procedure of integrating the defining dif- 
ferential equation over [R — e, R+ e] and taking the limit 
e \ 0. In the Newtonian case, this gives [defining p_ as 



This estimate is not angle averaged, but Mannarelli, Rajagopal, 
and Sharma's calculation has relatively large uncontrolled re- 
mainders, so we do not worry about the effects of angle averaging 
here. 
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the density immediately inside the star's surface and R 
to mean evaluation at R — e in the limit e \ 0] 



4nG 



p- 



3 
R 



(50) 



and in the GR case, we have (with G = 1) 



H' (R ) = H' old (R) 



4irh 



P-H (R), (51) 



where H' M (R) is computed using Eq. (f39|) . We 
thus make the replacement 3RF(R) — > [3 — 
4ttG p-Rj 'g(R)]RF '(R) in the expression for the New- 
tonian Green function [Eq. ([23"]) ]. and the replacement 
H' {R) -> H^ old (R)+A7rhp-H (R)/ct)'(R) in the GRcase 
[Eq. These changes in the boundary conditions in- 

crease the maximum quadrupole by a factor of < 2 in 
the example case considered below; the largest effect is 
for the least massive stars considered. 

Second, we would have to keep the boundary terms at 
the outer boundary when integrating by parts to obtain 



the expressions for the maximum quadrupole, since the 
shear modulus no longer vanishes at the star's surface. 
However, since here the shear modulus is smooth, it is 
numerically preferable not to perform any integration by 
parts, thus avoiding potential problems with large can- 
cellations between the surface and integrated terms. In 
this case, the expressions for the quadrupole assuming 
the UCB maximum uniform strain are [cf. Eqs. (|19[) and 
(131 



i^UCB strain, N i 
I W22 



32tt 
~l5~ 



and 



G N (r)[rp"(r) - p'(r)]dr 
(52) 



i^UCB strain, GR i 



where 



32tt 



15 J ^m^r) (53) 
+ G GR (r)4f B (r)R, 



T UCB / n 



V'-l)-2^ + ^ , + ^ 1 -^-^ 
r h 1 / 1 



p(r) 



2 + r(f - VQ 
h 1 ' 2 



r UCB 
'[t] 



(r) := 2 



p'(r) 



rp"{r) 
h 1 / 2 



b' + 2ib') - 5] - rib 1 - r 2 <b" - 1 , , , /? ) , s fr6' 
^ ^ 1 " + r0' + h 1 ' 2 \ p{r) - 2r f + 1 ) //(r). 



(54a) 
(54b) 



r 



However, these expressions will not actually yield the 
maximum quadrupole in this case, due to an important 
difference between the cases where the shear modulus 
vanishes at the star's surface and those where it does 
not. It is simplest to see this in the Newtonian case for a 
star with a constant shear modulus: Since the UCB max- 
imum strain expression (|52p only depends upon deriva- 
tives of the shear modulus, it predicts a zero maximum 
quadrupole, which seems absurd. One can, however, 
make a small adjustment to the form of the maximum 
strain one considers to yield a nonzero quadrupole in this 
case. This modification will also yield considerably larger 
maxima in the realistic case we consider, as well, where 
the shear modulus is close to constant — it decreases by 
less than a factor of 2 in going from the star's center to 
its surface in the example case we consider below. 

Specifically, in the case of a slowly varying shear mod- 
ulus, with /i(r) \rp'(r)\, \r 2 p"(r)\, appropriate for 
strange quark stars, we want the terms involving p, it- 
self to be largest. The appropriate choice for the strain 
in this case is most readily apparent from inspection of 
the Newtonian expression for the maximum quadrupole 
in terms of the stress tensor components, Eq. (|15p . We 
want the maximum contribution from the undifferenti- 
ated terms, which implies that we want t rr and — t r ± 



to be as large as possible. For <a, we note that since 
p'(r) < 0, we also want — to be as large as possible. 
Realizing that we can freely change the sign of any of the 
a. that give maximum uniform strain [given for the New- 
tonian case in Eqs. (fT"7)) ; cf. Eq. (65) in UCB], we thus 
reverse the sign of a r ± and a a ■ [The same logic holds for 
the more involved GR well, where the appropri- 

ate expression for a a will be the negative of Eq. (I3T)1) .] 

The resulting expressions for the putative maximum 
quadrupole in this case are thus 



i^-jmod. strain, I 
1^22 



^f/ GjV(r) [v Mr)+5/i ' (r) 

rp"(r) 



dr 



(55) 



and 



\Q 



mod. strain, GRi 
22 



-Vfr I [^m^m (56) 

+ Gc,E(.-)IH od (r)]*, 



where 
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j-mod 



[6p 7 5p]( r ) : 



M(r) + 



2 + r^' - V>') 



/i'(r) 



r/i"(r) 



r<f>'[5 - r(0' + 2ip')} + r^ 1 + r 2 <j>" + 1 



r<f>' + h 1 / 2 } l i(r)-2r + ^' (r) - 



In principle, these merely give a lower bound on the max- 
imum quadrupole, unlike the case in which the shear 
modulus vanishes below the surface, where there is a firm 
argument that maximum uniform strain maximizes the 
quadrupole. However, even if they do not give the abso- 
lute maximum, they should be quite close for cases like 
the one we consider here, where the shear modulus varies 
quite slowly. 
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FIG. 11: The Newtonian Cowling, Newtonian no Cowling, 
and full relativistic (including stress contributions) values for 
the quadrupole deformations (and fiducial ellipticity) of max- 
imally strained strange quark stars versus mass, using the 
EOS discussed in the text with a breaking strain of 0.1. We 
show these both for the standard UCB uniform maximum 
strain, and our modification that yields significantly larger 
quadrupoles in this case. 

Applying these expressions to a specific case, we use 
the strange quark matter EOS calculated by Kurkcla, Ro- 
matschke, and Vuorinen (KRV) [30| . generating an EOS 
for the parameter values of interest using the Mathe- 
matica notebooks available at The relevant pa- 

rameters are the values of the MS renormalization point, 
Aj^g and the strange quark mass, m s , both at a scale of 
2 GeV, along with the coefficient in the relation between 
the renormalization scale and the quark chemical poten- 
tial, X, the color superconductivity gap parameter, A 
(taken to be independent of density), 6 and the minimal 
quark chemical potential at which strange quark matter 



Note that A enters the KRV EOS through a color flavor locked 
(CFL) pressure term. This is not quite appropriate for the crys- 
talline color superconducting phase we consider here, since it 



(57a) 
(57b) 

I 

exists, /Ug, m in- We consider the EOS obtained by choosing 
A^ = 355 MeV, m s = 70 McV, X = 4, A = 10 McV, 
and fig^min — 280 McV. This parameter set yields a max- 
imum mass of 2.45 M@, with a maximum compactness of 
0.467. (These parameter choices were generally inspired 
by those considered at [(35[ , though with a smaller value of 
A, to place us well within the crystalline superconducting 
regime.) We show the quadrupole for a maximally uni- 
formly strained star in the three approximations (Newto- 
nian Cowling, Newtonian no Cowling, and GR) for both 
the UCB and modified maximum strain choices for this 
EOS in Fig. QTJ Here we have used a breaking strain of 
0.1, by the same high pressure argument as in the mixed 
phase case. (While the very outermost portions of the 
star are at low pressure, the parts that are at a lower 
pressure than the crustal case for which the 0.1 break- 
ing strain was calculated make negligible contributions 
to the quadrupole.) 



N.B.: To obtain the EOS used for this figure, we made 
some slight modifications to the KRV EoScalc Math- 
ematica notebook so that it would output particle num- 
ber densities on a denser mesh for low strange quark 
chemical potentials. This then gave an EOS table with 
better low-pressure coverage than their default settings 
produced. We still needed to perform an extrapolation 
of the EOS to zero pressure, where we found that a lin- 
ear extrapolation of the energy density and quark chem- 
ical potential in terms of the pressure using the lowest 
two entries of the table provided a good fit. (More in- 
volved approaches involving fitting to more points and/or 
a quadratic extrapolation produce very similar results.) 



assumes that all the quarks pair, while only some of them pair in 
the crystalline phase. However, as discussed in Sec. VI B of |66H , 
the condensation energy of the crystalline phases is easily 1/3 
to 1/2 that of the CFL phase with zero strange quark mass, 
which is the pressure contribution used by KRV. We have thus 
not altered this term in our calculations, since the contribution 
is already approximate, in that it assumes a density-independent 
gap parameter. Moreover, we only consider a fairly low value of 
A, while Knippel and Sedrakian |19| suggest that the crystalline 
phase might be favored up to A = 100 MeV. Our EOS may thus 
simply correspond to a slightly larger value of A, which would in- 
crease the maximum quadrupole, since the shear modulus scales 
as A 2 . 
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V. DISCUSSION 

Previous studies of the tidal and magnetic deforma- 
tions of compact stars have found similar rclativistic 
suppressions of quadrupole moments with compactness. 
In the tidal case, see the Love number computations 
in [H, [63|, HI, H3, HU- In the case of magnetic defor- 
mations, the expected suppressions are seen in, e.g., p~3l — 
[l5l |. In fact, since the largest compactness considered in 
these latter papers is only 0.48 (in [l|), one imagines 
that they overestimate the maximum quadrupoles by at 
least a factor of a few for more compact stars (for a fixed 
magnitude of magnetic field). 

As was argued by Damour and Nagar [64[ in the tidal 
case, all these suppressions are primarily related to the 
"no-hair" property of black holes: The largest relativistic 
suppression we find comes from the boundary conditions 
[through the Hq(R) and Hq(R) in the Green function's 
denominator — see Eq. (f4"2"j) ] , where one matches on to the 
external vacuum spacetime. For instance, for the SLy 
EOS's maximum compactness of 0.6, Ho(R) and Hq(R) 
are ~ 3.5 and ~ 6 times their Newtonian values [which 
can be obtained from the first term of Eq. (21) in Hin- 
derer [39| ]. In fact, these ratios go to infinity in the for- 
mal black hole limit, where the compactness approaches 
unity, as required by the no-hair property, and discussed 
by Damour and Nagar [64| (see their Sees. IV C and 
VII A, but note that their definition for the compactness 
is half of ours). This implies that the stiffness of spheri- 
cally symmetric curved vacuum spacetime suppresses the 
quadrupole. The quadrupole is also suppressed by a 
larger effective gravitational acceleration (given by </>'), 
which appears in the denominator of Gqr, replacing the 
Newtonian g(r) [cf. Eqs. flU} and (|15|)]. (But recall that 
we always compute the background stellar structure rel- 
ativistically, so this larger acceleration only affects the 
perturbation equations, and not, e.g., the thickness of 
the crust for a given mass and EOS, which is the same in 
both the Newtonian and relativistic calculations of the 
quadrupole.) 

Our results imply that nearly all of the Newtonian 
computations of quadrupoles due to elastic deformations 
of relativistic stars overestimate the quadrupole moment, 
often by at least a factor of a few. The only exceptions 
we have found are for low-to-mid mass strange quark 
stars and for elastic stresses in the cores of neutron stars 
around 0.5 M Q . In both of these cases, the Newtonian 
Cowling approximation is a slight underestimate for con- 
tributions to the quadrupole, though the Newtonian no 
Cowling version is still an overestimate. See Fig. [T] for an 
illustration in the core case; but note that neutron stars 
with such low masses are not known to exist in nature. 
The overestimate from performing a Newtonian Cowl- 
ing approximation calculation can be ~ 6 for massive 
stars whose quadrupole is being generated by an elastic 
deformation near the crust-core interface, as considered 
by UCB and others. This is due in part to the sud- 
den changes in density at that interface entering directly 



through g", as discussed at the end of Sec. IPTl 

However, the calculations by Horowitz [2(| for crustal 
deformations of very low mass stars only receive negli- 
gible corrections (of < 5%), since he considers compact- 
nesses of ~ 0.01. In fact, one makes even smaller errors 
in using the Cowling approximation to treat these stars, 
since the changes in density in the crust (times AnGr 2 ) 
are much smaller than the star's gravitational field there. 

No neutron stars with such low masses have ever been 
observed (nor is there a compelling mechanism for form- 
ing them). Nevertheless, Horowitz remarks that gravi- 
tational wave detection of gravitational waves from elas- 
tically deformed neutron stars will, ceteris paribus, be 
biased towards low(er) mass neutron stars, if one consid- 
ers deformations generated by crustal stresses. This is an 
important point, particularly when considering the astro- 
nomical interpretation of detections (or even upper lim- 
its), and the results we present here make the bias against 
high-mass stars even stronger. (This bias also applies to 
solid quark stars, though there it is rather weak. It does 
not apply to hybrid stars, however, where it is high-mass 
stars that can sustain the largest quadrupoles.) 

Of course, one must remember that all of these val- 
ues are maxima, assuming a maximally strained star, 
while there is no reason, a priori, for a given star to be 
maximally strained. Moreover, as UCB and HJA note, 
these calculations assume that all the strain goes into the 
I = m = 2 perturbation, though strain in other modes 
(e.g., the I = 2, m = mode due to rotation) can push 
the lattice closer to its breaking strain while not increas- 
ing the I = m = 2 quadrupole. 



VI. CONCLUSIONS AND OUTLOOK 

We have presented a method for calculating the maxi- 
mum elastic quadrupole deformation of a relativistic star 
with a known shear modulus and breaking strain. We 
then applied this method to stars whose elastic deforma- 
tions are supported by a shear modulus cither from the 
Coulomb lattice of nuclei in the crust, a hadron-quark 
mixed phase in the core, or crystalline superconducting 
strange quark matter throughout the star. (In the last 
case, we have made the requisite changes to the method 
so that it is valid when the star has a nonzero surface 
density and the shear modulus does not vanish at the 
star's surface.) In all but the strange quark case, we 
find that the relativistic quadrupole is suppressed, com- 
pared with the standard, Newtonian Cowling approxi- 
mation calculation of the quadrupole, at least for stars 
with masses of > 1 Mq (corresponding to the observed 
masses of neutron stars) and the EOSs we have investi- 
gated. These suppressions can be up to ~ 4 in the hybrid 
case, and ~ 6 in the crustal case. In the strange quark 
star case, the Newtonian Cowling approximation calcu- 
lation slightly underestimates the quadrupole (by tens of 
percent) for low-to-standard mass stars, but is still an 
overestimate of ~ 2 at higher masses. 
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These suppressions strengthen the Horowitz [2Cj ar- 
gument that searches for gravitational waves from elas- 
tically deformed neutron stars supported by crustal 
stresses are biased towards lower-mass stars. The same 
argument also applies to strange quark stars, though 
there the suppressions with increasing mass are less se- 
vere (and the maximum quadrupoles are all consider- 
ably larger). However, this argument does not apply to 
quadrupolc deformations of hybrid stars, since the in- 
crease in the size of the region of mixed phase with in- 
creasing mass dominates the various suppressions. 

Our results also imply that many of the pre vious cal- 
culations of elastic quadrupoles (e.g., fl7H19l . l2p . HU) 
will need their results revised downwards. (While we 
find much larger maximum quadrupoles for solid strange 
quark stars than did Lin [18j . this is only because we 
assume a breaking strain 10 times that assumed by Lin. 
If we take the same 10~ 2 breaking strain as does Lin, 
then we find a suppression of a factor of a few, though 
this is very likely within the uncertainties of Lin's calcu- 
lation, which assumed a uniform density, incompressible 
star with a uniform shear modulus.) 

It is instructive to compare our results with the num- 
bers quoted in Pitkin's review 0. All of these were 
obtained by Pitkin using scalings given in the aforemen- 
tioned pap ers, sometimes updating to the Horowitz and 
Kadau [3l| breaking strain, and provide a good overview 
of the standard Newtonian predictions. None of our de- 
tailed calculations for maximum crustal quadrupoles ap- 
proach the high values Pitkin obtained using UCB's fit- 
ting formula (as corrected by Owen [HI). However, our 
very largest hybrid star quadrupoles are an order of mag- 
nitude above Pitkin's quoted maximum, even if one only 
assumes a breaking strain of 10 -2 , as does Pitkin. Ad- 
ditionally, our estimates for maximum solid quark star 
quadrupoles (~ 10 44 g cm 2 for 1.4 Mq stars) are consid- 
erably larger than the ones quoted by Pitkin (based on a 
different shear modulus model), even if we reduce them 
by an order of magnitude due to scaling the breaking 
strain to Pitkin's 10~ 2 . In fact, they are in the same 
range as those Pitkin quotes for a model for crystalline 
superconducting hybrid stars (with an optimistic gap pa- 
rameter 5 times the one we used for solid quark stars, 
leading to a shear modulus ~ 40 times our shear modu- 
lus's maximum value). 

Even with the relativistic suppressions, we obtain max- 
imum quadrupole deformations of a few x 10 42 g cm 2 
in the hybrid case for a very stiff hadronic EOS, and 
a few x 10 41 g cm 2 for more realistic cases. In both sit- 
uations, the largest maximum quadrupoles are given by 
the most massive stars. These values are proportional 
to the brea king strain and assume that the Horowitz 
and Kadau [3l| breaking strain of about 0.1 is appli- 
cable to the mixed phase. Such large quadrupole defor- 
mations were previously thought only to be possible for 
solid quark stars (see [lH EH, EH H3 ) , or from crustal de- 
formations in the very low-mass neutron stars considered 
by Horowitz [20j . These large deformations (correspond- 



ing to fiducial cllipticities of a few x 10 -3 in the extreme 
case, and ~ 5 x 10~ 4 in a more realistic case) would be 
able to be detected by current LIGO searches for grav- 
itational waves from certain known neutron stars |6j-l8j. 
(However, we must note that there is no reason to as- 
sume that such isolated stars are anywhere near maxi- 
mally strained, even neglecting the uncertainties in the 
description of their interiors.) 

The prospects for crustal quadrupoles are now some- 
what less optimistic, and definitely favor lower-mass 
stars. However, for a canonical 1.4 Mq neutron star, we 
find that the maximum relativistic crustal quadrupole 
is in the range ~ (1-6) x 10 39 g cm 2 [corresponding to 
fiducial cllipticities of ~ (1-8) x 10 -6 ], depending on the 
model used for the crust and the high-density EOS. (Note 
that the fully consistent Douchin and Haensel model with 
its associated high-density EOS yields the lowest num- 
bers. Additionally, there is the possibility of a further 
reduction of up to ~ 2 due to the angle averaging pro- 
cedure used to obtain the shear modulus.) On the high 
side, these numbers are consistent with those given pre- 
viously for breaking strains of 0.1 by Horowitz [2(| l3l| . 
though they are a factor of ~ 5 lower than the maximum 
Pitkin [l(| obtained using scalings of previous results and 
the maximum value given by H J A (scaled to this breaking 
strain). For stars around 2M Q , the relativistic suppres- 
sions lead to maximum quadrupoles that are nearly an or- 
der of magnitude smaller than those for a 1.4 Mq star in 
the compact SLy case: ~ (1-5) x 10 38 g cm 2 [correspond- 
ing to fiducial ellipticitics of ~ (1-6) x 10~ 7 ]; and even in 
the much less compact LKR1 case, there is a suppression 
of ~ 5. Previous Newtonian studies (see Fig. 3 in [2(|) 
had only found suppressions of around a factor of 4, due 
to the thinning of the crust and the increase in Newtonian 
gravity with increasing mass. It will be interesting to 
consider further models for the crustal composition and 
EOS in this case, particularly the large suite of crustal 
models including the pasta phases recently calculated by 
Newton, Gearheart, and Li j6{|. (See [7fl| for order- 
of-magnitude estimates of the maximum quadrupole for 
these models, illustrating the sensitive dependence on the 
slope of the symmetry energy.) 

One would also like to make similar relativistic calcula- 
tions of the maximum energy stored in these elastic defor- 
mations. This would be useful in properly computing the 
available energy for magnetar flares, for instance. (Using 
Newtonian scalings, Corsi and Owen [7l[ estimated that 
the hybrid case was especially interesting compared to 
existing LIGO upper limits for gravitational wave emis- 
sion from such flares.) The basic expressions (at least 
in the perfect fluid case) a ppe ar to be readily available 
in the literature (see, e.g., |72l. [73j ; [H, EH give related 
results including elasticity). However, one cannot apply 
these directly to the crustal and hybrid cases, even in 
the Newtonian limit, due to the distributional nature of 
the density and pressure perturbations. Specifically, the 
sudden change in shear modulus at the phase transitions 
gives delta functions in the derivatives of the density and 
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pressure perturbations. Since the energy expressions in- 
volve squares of these derivatives, one would have to in- 
voke some sort of rcgularization procedure, or apply a 
different method. Developing appropriate expressions for 
this case will be the subject of future work. 

Returning to the quadrupoles, one might also want 
to consider the shape of the deformed star, particularly 
in the relativistic case — the ellipticity is already only a 
rough indicator of the shape of the deformation in the 
Newtonian case. But the effects of the star's magnetic 
field are surely the most interesting to consider, from 
its influence on the lattices that support elastic defor- 
mations, to the changes to the boundary conditions at 
the star's surface from an external magnetic field (par- 
ticularly for magnetars), to the internal magnetic field's 
own contribution to the star's deformation. One might 
also want to consider the lattice's full elastic modulus 
tensor in this case, instead of simply assuming a poly- 
crystalline structure and angle averaging to obtain an ef- 
fective isotropic shear modulus, as was done here. (And 
even if one assumes a polycrystalline structure, one could 
use more involved, sharper bounds on the shear modu- 
lus than the ones considered here — see (2fl for a classic 
review of such bounds.) 
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Appendix: Hadron— quark hybrid stars and the 
binding energy argument 

As we mentioned in Sec. II C of [28(, if the surface 
tension is large enough, the mixed phase is not locally 
favored energetically (i.e., at a fixed baryon density), 
compared to the individual pure phases. (This was first 
noted by Heiselberg, Pethick, and Staubo [76| and later 
discussed by Alford et al. [23].) However, as was also 
noted in [281 ] . the entire region of mixed phase can still 
be favored due to global energy arguments, especially 
when one considers the binding energy of the star (for a 
fixed total baryon number): One expects the stars with 
the largest binding energy (i.e., smallest gravitational 
mass) for a given total baryon number to be favored. 
Here we always compare with a purely hadronic star. 
One would expect the Maxwell construction case with a 
sharp interface between the two phases to produce more 
strongly bound stars, given the local energy results pre- 
sented in [76|,[23]. However, at least for the EOSs consid- 
ered in Fig. [T21 the Maxwell construction stars with total 
baryon numbers up to that of the corresponding maxi- 
mum mass hybrid star only contain hadronic matter. 




-0.001 - 
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FIG. 12: The differences in hybrid and hadronic star masses 
versus mass for a fixed total baryon number versus hybrid 
star mass for a variety of equations of state and surface ten- 
sions. These calculations include the lattice contributions to 
the EOS except in the "plain" cases. 



Specifically, if we compute the gravitational mass of a 
star with a given total baryon number using, e.g., the 
Hyl EOS parameters, we find that this mass is smaller 
(corresponding to a larger binding energy) when includ- 
ing the mixed phase than when considering a purely 
hadronic star (constructed with the same hadronic EOS 
parameters). However, the differences in these masses 
are not very large, and the situation is reversed for lower 
masses — and often for all masses — if one includes the con- 
tributions of the lattice to the EOS, viz., the lattice's en- 
ergy density and pressure, and the contributions of the 
surface tension to the energy density. [We compute the 
first two contributions using Eqs. (14) and (20) in [28| . 
and the last contribution using that paper's Eq. (2).] 
This is illustrated for a few representative EOSs and sur- 
face tensions in Fig. [T^J some relevant quantities for stars 
constructed with these parameters are given in Table [U 
(We calculate these mass differences by first computing 
the total baryon number as a function of mass for the 
purely hadronic stars and then using bisection to locate 
the hybrid star with the same total baryon number.) 

The resulting mass differences severely restrict the 
range of masses for which the mixed phase exists for 
reasonable EOS parameters, and often eliminate it com- 
pletely, at least for the EOSs we have considered. Indeed, 
of the EOSs considered in [28| (and for surface tensions 
for which all of the lattice is stable), the mixed phase 
is only favored in stable stars for the two more extreme 
cases, viz., Hyl' and LKR1, and there only for low sur- 
face tensions, as shown in Fig. [T2] Additionally, in these 
cases, the binding energy argument would imply that one 
would suddenly nucleate a large region of mixed phase if 
one increased the star's mass to just above the mass for 
which the mixed phase is first favored — this sort of pro- 
cess is discussed in Chap. 7 of [48|. 

However, as was also mentioned in [28| . our energy 
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a 


EOS additions M max 


■ .-hybrid 
mill 


phybrid 
-'••max 


/ R Cmax 


densest 




MeV far 2 






Mq 


% 




hybrid phase 






none 


1.955 


1.10 


72.4 


0.433 


H, d = 3.00 


LKR1 




blob + lattice 


1.944 


1.04 


( 1.6 


n a q i 
U.4ol 


tt j o nn 
tl, a = o.UU 


20 


blob + lattice 


1.937 


1.79 


72.3 


0.430 


H, d = 3.00 




30 


blob + lattice 


1.930 


1.87 


72.2 


0.430 


H, d = 3.00 


Hyl' 




none 


1.973 


1.38 


69.3 


0.476 


H, d= 1.51 


20 


blob + lattice 


1.961 


1.89 


68.8 


0.474 


H, d= 1.30 



TABLE I: Extent of the mixed phase in stable stars for two EOSs from (2^], with various surface tensions. In the "EOS 
additions" column, "blob + lattice" denotes the case where we have included the blob energy density and lattice energy density 
and pressure in the EOS. The other columns are the same as in Table I of H|. Explicitly, Af££ rid gives the masses of stars that 
first contain hybrid matter (using the binding energy argument); Rm^ 1 /R denotes the maximum radius fraction occupied by 
hybrid matter (i.e., the radius fraction for the maximum mass star); and C max denotes the maximum compactness (2GM/Rc 2 ) 
of a star. We also give the composition of the rare phase ("H" stands for hadronic) and the dimension of the lattice at the 
center of the maximum mass star. 



computations are sufficiently uncertain that these bind- 
ing energy results should not be taken as definitively rul- 
ing out the mixed phase for these parameters, since the 
differences in binding energies (and EOSs) in question 
are not large. In particular, a self-consistent treatment of 
charge screening would likely result in lower lattice ener- 
gies, and thus make the mixed phase favored for a larger 
range of par ameters. Indeed, Christiansen and Glenden- 
ning [78l. l79j| argue that the mixed phase should always 
be favored, and any calculation that predicts otherwise 
must be incomplete or using inapplicable input parame- 
ters. It thus seems reasonable to treat the plain results 
presented in Fig. |H] and those including the lattice con- 
tributions to the EOS and using the binding energy ar- 
gument as roughly bracketing the possible values for the 
maximum quadrupole. (However, the binding energy ar- 
gument results are not a firm lower bound.) 

N.B.: We have not shown the calculation for lower sur- 
face tensions since the lower-dimensional portions of the 
lattices are unstable for surface tensions not too much 
smaller than the ones we consider here. In particular, the 



lower-dimensional portions of the lattice are unstable for 
the Hyl and Hyl' EOSs for a = 10 MeV fm~ 2 and for 
the LKR1 EOS for a = 5 MeV fm~ 2 . All this may mean 
is that there is a different lattice geometry that would 
actually occur for those surface tensions, but we can- 
not make predictions for such cases with the calculations 
from [28| , so we simply exclude them from consideration. 

We have illustrated the quadrupoles for a variety of 
lower surface tensions in Fig. 1101 where we show the 
minimum masses for which hybrid matter is present ac- 
cording to the binding energy argument. While we have 
not shown this in the plot, including the lattice contri- 
butions to the EOS increases the maximum quadrupole 
by ^ 25% for the surface tensions we have considered. 
However, for large surface tensions (e.g., 80 MeV fin ), 
including these contributions also reduces the maximum 
mass to the point where it is no longer compatible with 
the Demorest et al. [8(| measurement of a 1.97±0.04M Q 
neutron star for EOSs which already have a rather low 
maximum mass (e.g., LKR1 and Hyl'). 
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